Radiation-induced magnetoresistance oscillations in two-dimensional electron systems 

under bichromatic irradiation 
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We analyze the magnetoresistance R xx oscillations in high-mobility two-dimensional electron sys- 
tems induced by the combined driving of two radiation fields of frequency uj\ and 0J2, based on 
the balance-equation approach to magnetotransport for high-carrier-density systems in Faraday ge- 
ometry. It is shown that under bichromatic irradiation of U2 ~ 1.5u>i, most of the characterstic 
peak-valley pairs in the curve of R xx versus magnetic field in the case of monochromatic irradiation 
of either uj 1 or 102 disappear, except the one around wi/w c ~ 2 or wa/Uj ~ 3. R xx oscillations show 
up mainly as new peak-valley structures around other positions related to multiple photon processes 
of mixing frequencies u> 1 + u>2, ^2 — Wi, etc. Many minima of these resistance peak- valley pairs can 
descend down to negative with enhancing radiation strength, indicating the possible bichromatic 
zero-resistance states. 
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PACS numbers: 73.50.Jt, 73.40.-c, 78.67.-n, 78.20.Ls 



I. INTRODUCTION 



Since the discovery of radiation induced mag- 
netorersistance oscillations (RIMOs) and zero- 
resistance states (ZRS) in ultra-high mobil- 
ity two-dimensional (2D) electron systems 
tremendous experimenta l 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 ' 13 ' 14 and 
theoretica l 15 ' 16 - 17 - 18 - 19 - 20 - 21 - 22 - 23 - 24 - 25 - 26 - 27 - 28 - 29 - 30 - 31 - 32 - 33 - 34 
efforts have been devoted to study this exciting phe- 
nomenon and a general understanding of it has been 
reached. Under the influence of a microwave radiation of 
frequency / = uj/2tt, the low-temperature magnetoresis- 
tance R xx of a 2D electron gas (EG), exhibits periodic 
oscillation as a function of the inverse magnetic field. 
The RIMOs feature the periodical appearance of peak- 
valley pairs around lo/lj c = 1, 2, 3, 4, • • •, i.e. a maximum 
at uj/u) c = j — 5 J and a minimum at lu/lu c = j + 5j, 
with j = 1, 2, 3, 4, • • • and < Sf < 1/4. Here uj c is the 
cyclotron frequency and lu/lo c = j are the node points of 
the oscillation. With increasing the radiation intensity 
the amplitudes of the peak-valley oscillations increase 
and the resistance R xx around the minima of a few 
lowest-j pairs can drop down towards negative direction 
but will stop when a vanishing resistance is reached, i.e. 
ZRS. 

In addition to these basic features, sec- 
ondary peak-valley pair structures were also ob- 
served experiment ally^&LSiiiii^ and predicted 
theoretically^! around uj/uj c = 1/2,3/2 and 2/3. 
They were referred to the effect of two- and three-photon 
processes and their minima were shown also to be able 
to develop into negative value when increasing the 
radiation further. 23 Recent measurement at 27 GHz 
with intensified microwave intensity^ confirmed these 
mutliphoton-related peak-valley pairs and the ZRSs 
developed from their minima. A careful theoretical anal- 
ysis with enhanced radiation convincingly reproduced 
these structures and predicted more peak-valley pairs 
related to multiphoton processes^ 



Despite intensive experimental and theoretical stud- 
ies have been done in the case of monochromatic irra- 
diation, further investigations beyond this configuration 
are highly desirable for a deeper understanding of this 
fascinating phenomenon. An easy way is to study the 
response of the 2D system to a bichromatic radiation, 
which apparently can not be reduced simply to the super- 
position of the system response to each monochromatic 
radiation. Thus the presence of a second radiation of dif- 
ferent frequency could provide additional insight into the 
problem of microwave-driven 2D electron system. 

Theoretically, effect of a bichromatic irradiation on 
transport of a 2D electron gas was investigated within 
a model of a clean classical gas, in which Kohn's theo- 
rem is violated entirely due to nonparabolicity without 
invoking Landau quantization. 36 Some features that are 
specific to the bichromatic case, such as new domain of 
magnetic field within which the diagonal conductivity is 
negative, were revealed. 

Experimentally, Zudov et alH recently measured the 
magnetoresistance of a high-mobility 2D electron system 
under the combined driving of two radiation fields of fre- 
quency u>i and u>2, and disclosed features of RIMOs which 
are quite different from those under monochromatic ra- 
diation of frequency lj± or L02. They detected a new resis- 
tance minimum under bichromatic microwave radiation, 
which seems to originate from a frequency mixing pro- 
cess, possibly a precursor of bichromatic ZRS. 

In this paper we report our studies on microwave pho- 
toresistance response of high-mobility two-dimensional 
electron systems under bichromatic irradiation, based 
on the balance-equation approach to magnetotransport 
for high-carrier-density systems, extended to the case 
of simultaneous driving of two radiation fields of differ- 
ent frequencies. The balance-equation approach, though 
semiclassical in nature, has been shown to capture the 
essence of this radiation-induced nonlinear magneto- 
transport quantitatively^^! Under bichromatic irradi- 
ation of u>2 ~ 1.5wi, we find that most of the characteris- 
tic peak- valley pairs in the curve of R xx versus magnetic 
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field in th case of monochromatic irradiation of either 
u)% or L02 disappear, except the one around u)\/u) c ~ 2 
or w 2 /cj c ~ 3. i?;^ oscillations show up mainly as new 
peak-valley structures around other positions related to 
multiple photon processes of mixing frequencies u>i + cj 2 , 
wa — £Ji, etc. 

II. BALANCE EQUATIONS UNDER 
BICHROMATIC RADIATION 

The derivation of balance equations under bichromatic 
radiation follows that of monochromatic radiation. 31 The 
nature of the balance-equation approach and its appli- 
cability to radiation-driven magnetotransport in high- 
mobilty 2D electron systems of high carrier density, was 
examined closely in Ref.^l We refer the readers to it for 
detail. 

In the case of bichromatic radiation we consider that a 
dc or slowly time- varying electric field Eo and two high 
frequency (HF) fields 

Ei(£) = Ei s sin(wii) + E ic cos(cJii) (1) 

and 

E 2 (t) = E 2s sin(w 2 t) + E 2c cos(w 2 i) (2) 

are applied simultaneously in a quasi-2D system consist- 
ing of iV c interacting electrons in a unit area of the x- 
y plane, together with a magnetic field B = (0, 0, B) 
along the z direction. The frequencies lu± and a-> 2 are 
high enough and their difference is large enough that oj\ 
and u>2, as well as — n 2 cj 2 | (for arbitrary integers n\ 
and n 2 ), are all much larger than 1/tq, where tq stands 
for the scale of the time variation of slowly varying field 
Eo, or the time scale within which one carries out the 
transport measurement, whichever the shorter. The ap- 
proach is based on the separation of the center-of-mass 
motion from the relative electron motion of the electrons 
and describes the transport state of a high-carrier-density 
many-electron system under radiation fields in terms of 
a rapidly time-varying electron drift velocity oscillating 
at both base radiation frequencies, v(i) = vi(i) + v 2 (t), 
with 

Vi(i) = vie cos(wit) + Vis sm(wit), (3) 

v 2 (i) = v 2c cos(u> 2 i) + v 2s sin(a; 2 i), (4) 

together with another part Vo for the slowly varying elec- 
tron drift motion, as well as an electron temperature T e 
characterizing the electron heating. 37,38 In the case of 
ultra-clean electron gas at low temperatures, the slowly 
time-varying quantities Vo and T e satisfy the following 
force- and energy-balance equationsiSii^i 

m -jT = eE ° + e ( v n x B ) + tt> ( 5 ) 



iV c E • v + S p - W = 0, (6) 

with Vic and vi s determined by 

- mwivic = eEi s + e(vi s x B), (7) 

mwivi s = eEi c + e(v ic x B); (8) 

and v 2c and v 2s determined by 

-mw 2 v 2c = eE 2s + e(v 2s x B), (9) 

mw 2 v 2s = eE 2c + e(v 2c x B). (10) 



Here e and m are the electron charge and effective mass, 

00 

F = ^|c/( q|l )| 2 J2 (60 

q|| rii,n 2 = — 00 

x q|| II 2 (q|| , lo - mui - n 2 w 2 ) (11) 

is the damping force of the moving center-of-mass, 

00 

s P = £|£%n)l 2 E J 2 m&)Jl(& 

qn ni,na=— 00 

x (niwi + n 2 LU 2 ) n 2 (q||,wo - - «2^ 2 ) (12) 

is the averaged rate of the electron energy absorption 
from the HF fields, and 

00 

VF = ]T|M(q)| 2 ]T ^(£1)^(6) 

q ni,7i2 — — 00 

x O q A 2 (q, cj + f2 q -niwi — n 2 w 2 ) (13) 

is the average rate of the electron energy dissi- 
pation to the lattice. In the above equations, 
Jni(£i) an d J n2 (&) are Bessel functions of order m 
a nd n 2 , £1 = ^/fqjj ■ Vi e ) 2 + (q|| ■ w 1s ) 2 /uji, £> = 

\/(q|| • v 2c ) 2 + (q|| • v 2s ) 2 /^2, and u = qy • v . Here 
9|| = {Qx, Qy) stands for the in-plane wavevector, f7(qii) is 
the effective impurity scattering potential, and n 2 (q| , O.) 
is the imaginary part of the electron density correlation 
function of the quasi-2D system in the magnetic field. 
In Eq. (|13f) . q represents the 3D wavevector (qy , q z ) plus 
the branch index A, and the summation is for all possi- 
ble 3D phonon modes of frequency Qq having electron- 
phonon scattering matrix element A/(q). A 2 (q, ft) = 
2n 2 (q||,n)[n(n q /T)-n(n/T e )] (with n(x) = l/(e x -l)) 
is the imaginary part of the electron-phonon correla- 
tion function. The n 2 (q||,f2) function of a 2D sys- 
tem in a magnetic field can be expressed in the Landau 
representation* 3 ^ 

n 2 (q| h ^) = ^2^C„, n ,(/^2/2)n 2 (n,<0), (14) 

n 2 (n, n',!l)=~J de [/(e) - f(e + fi)] 

xImG„(£ + fi)ImG n /(e), (15) 
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where C n>n+ i(Y) = n![(n + l)^Y l e^ Y [L l n {Y)] 2 with 
L L n (Y) the associate Laguerre polynomial, Ib = \Jl/\eB\ 
is the magnetic length, f(e) — {exp[(e — u)/T c ] + 
is the Fermi distribution function at electron tempera- 
ture T e . The density of states of the n-th Landau level is 
modeled with a Gaussian form42> 

ImG„(e) = -(VVr)exp[-2(e - e n f /T 2 ], (16) 

having a half- width 

around the level center e n = tilu c . Here lu c — eB/m is 
the cyclotron frequency, /io is the linear mobility at lat- 
tice temperature T in the absence of the magnetic field, 
and a is a semi-empirical parameter to take account of 
the difference of the transport scattering time r m deter- 
mining the mobility uo, from the single particle lifetime 
t s related to Landau level broadening. 

For time-independent Vo, we immediately deduce the 
transverse and longitudinal dc resistivities from Eq. 0: 

R xy = B/N c e, (18) 



under given incident radiation fields, before directly cal- 
culating the linear magnetoresistivity from Eq. (|20[l . 

Within certain field range, the magnetoresistivity R xx 
given by Eq. (|19|1 can be negative at small i>o, but will 
change towards the positive direction with increasing vq 
and passes through zero at a finite v o implying that the 
time-independent small-current solution is unstable and 
a spatially nonuniform^ or a time-dependent solutionis 
may develop, which exhibits measured zero resistance. 
Therefore we identify the region where a negative dissi- 
pative magnetoresistance develops as that of the ZRS. 

The summations over n\ and 712 in Eqs. (|llfl . i|12|) and 
(fLl|l include all possible electron transition processes as- 
sisted by real and virtual photons of frequency ui\ and 
L02. The ri\ > terms represent all possible electron 
transitions with simultaneous absorption of n\ photons 
of frequency oj\, while ri\ < terms represent all possi- 
ble electron transitions with simultaneous emission of \n\ \ 
photons of frequency uj\. We call these electron transi- 
tions the I ni| uj\ -photon assisted processes. The same for 
|^2| W2-photon assisted processes. The m = terms rep- 
resent electron transitions assisted by virtual wi-photons 
(all possible emission and absorption of same number of 
cji-photons), and n-i = terms represent electron transi- 
tions assisted by virtual <x>2-photons. 



R xx = -F ■ v /{N 2 e 2 v 2 Q ). 
The (linear) magnetoresistivity is its vo — > limit: 

<9n 2 



(19) 



ri\ .712 — — 00 



on 



(20) 



We assume that the 2DEG is contained in a thin sam- 
ple suspended in vacuum at plane z = 0. When both 
electromagnetic waves illuminate upon the plane per- 
pendicularly with the incident electric fields Eji(t) = 
Eii s sin(wit) + Eii c cos(cjit) and Ei2(£) = Ej2 S s\n(u)2t) + 
E;2c cos(u>2t) 1 the HF electric fields in the 2DEG, deter- 
mined by the electrodynamic equations, are 

„ . . N c evi(t) , . . . 

Ei* = ° — +E il (t), 21 
2e c 

E 2 (*) = ^) + E i2 (t). (22) 

Using this E x (t) in Eqs. Q and ©, and E 2 (i) in Eqs. @ 
and IjlOp. the oscillating velocity vi c and Vi s (and thus 
the argument £1) are explicitly expressed in terms of 
incident field E;i s and Eji c , and the oscillating veloc- 
ity V2 C and V2 S (and thus the argument £2) are explic- 
itly expressed in terms of incident field Ej2 S and Ej2 C . 
Therefore, in the case of weak measuring current limit 
(vq — » 0) we need only to solve the energy balance equa- 
tion S„ — W = to obtain the electron temperature T e 



III. NUMERICAL RESULTS 

As indicated by experiments^ in ultra-clean GaAs- 
based 2D samples having mobility of order of 10 3 m 2 /Vs, 
the remote donor scattering is responsible for merely 
~ 10% or less of the total momentum scattering rate. 
The dominant contribution to the momentum scatter- 
ing comes from short-range scatterers such as residual 
impurities or defects in the background. Therefore, as 
in Ref.l3ll we assume that the dominant elastic scatter- 
ings contributing to the resistance and energy absorption 
are due to short-range impurities randomly distributed 
throughout the GaAs region in the numerical calcula- 
tions. The impurity densities are determined by the 
requirement that electron total linear mobility at zero 
magnetic field equals the giving value at lattice temper- 
ature T. Since RIMO measurements are at low temper- 
atures, the direct phonon contributions to S p and Fo 
(Rxx) can be neglected. Nevertheless, to calculate the 
electron energy dissipation to the lattice, W, we take ac- 
count of scatterings from bulk longitudinal acoustic (LA) 
and transverse acoustic (TA) phonons (via the defor- 
mation potential and piezoelectric couplings), as well as 
from longitudinal optical (LO) phonons (via the Frohlich 
coupling) in the GaAs-based system. The relevant mate- 
rial and coupling parameters are taken typical values of 
GaAsi^ electron effective mass m = 0.068 m e (m e is the 
free electron mass), transverse sound speed v s % = 2.48 x 
10 3 m/s, longitudinal sound speed v s \ = 5.29 x 10 3 m/s, 
acoustic deformation potential 5 = 8.5 eV, piezoelec- 
tric constant ei4 = 1.41 x 10 9 V/m, dielectric constant 
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FIG. 1: (Color online) The magnetoresistivity R xx versus 
the inverse magnetic field 1/B for a GaAs-based 2DEG with 
iV e = 3.0 x 10 15 m" 2 , fj, = 2000m 2 /Vs and a = 5, irradi- 
ated simultaneously by to two microwaves with frequencies 
u)i/2n = 31 GHz and uj2/2tt = 47 GHz having four sets of in- 
cident amplitudes En = Ei2 — 1.4, 2, 3, and 4 V/cm at lattice 
temperature T — 1 K. uo c = eB/m is the cyclotron frequency 
at the magnetic field B. 
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FIG. 2: (Color online) The magnetoresistivity R xx versus the 
inverse magnetic field 1/B for the same system as described in 
Fig. 1. The thin curves are under monochromatic irradiation 
either of frequency u)\/2-k — 31 GHz and incident amplitude 
Ei = 3 V/cm (dot curve), or of frequency u>i/2ir — 47 GHz 
and incident amplitude E\ = 3 V/cm (solid curve). The 
thick curves are under bichromatic irradiation of frequencies 
U3i/2n — 31 GHz and 0J2/2tt = 47 GHz having incident ampli- 
tudes En — Ei2 = 3 V/cm (solid curve) or 3.5 V/cm (dash-dot 
curve). The lattice temperature is T = 1 K. 



k = 12.9, material mass density d = 5.31 g/cm 3 . 

The numerical calculations are performed for in- 
direction (parallel to Eo) linearly polarized incident mi- 
crowave fields [Eii s = (En, 0), E;i c = and Ei 2;5 = 
OEk,0),Eac = 0]. 

Figure 1 shows the calculated magnetoresistivity R xx 
versus the inverse magnetic field for a GaAs-based 2D 



co,/2n=31 GHz (0/271=47 GHz 
T=1 K 




FIG. 3: (Color online) The same as Fig. 2, but magnetoresis- 
tivity R xx is shown as a function of the magnetic field B. 



system having electron density N c = 3.0 x 10 m , lin- 
ear mobility fig = 2000m 2 /Vs and broadening parameter 
a = 5, simultaneously irradiated by two microwaves of 
frequencies uji/2tt = 31 GHz and lo^I^ = 47 GHz with 
four sets of incident amplitudes En = E i2 = 1.4,2,3, 
and 4 V/cm at lattice temperature T = 1 K. We see that 
the bichromatic photoresponse of R xx is unlike that of 
monochromatic irradiation. Most of the characterizing 
peak-valley pairs in the curve of R xx versus magnetic 
field subjected to monochromatic irradiation of either wi 
or u)2 disappear. The only exception is the peak-valley 
pair around w\/u) c ~ 2 or u^s/lo c ~ 3, which, on the con- 
trary, is somewhat enhanced in bichromatic irradiation. 
New peak-valley pairs appear. Their amplitudes gener- 
ally increase with increasing the incident field strengths 
as indicated in the figure. Most of these pairs exhibit es- 
sentially fixed node positions while growing amplitudes 
and their minima can drop down into negative. 

As in the case of monochromatic illumination, the ap- 
pearance of magnetoresistance oscillation in the bichro- 
matic irradiation comes from real photon-assisted elec- 
tron transitions between different Landau levels as in- 
dicated in the summation of the electron density- 
correlation function in Eq. Ijl4jl. Apparently, all multi- 
ple photon processes related to w\ monochromatic radi- 
ation and LU2 monochromatic radiation are included. In 
addition, there are multiple photon processes related to 
mixing u>\ and L02 radiation. 

With three positive integers 71 1 , TI2 and I to characterize 
a multiple photon assisted electron transition, we use the 
symbol ri\U!\ + 712^2 : I to denote a process during which 
an electron jumps across / Landau level spacings with 
simultaneous absorption of 711 photons of frequency ui\ 
and n 2 photons of frequency 0J2, or simultaneous emis- 
sion of n% photons of frequency u>x and ri2 photons of 
frequency L02 ; and use the symbol ri\uj\ — 712W2 : I to 
denote a process during which an electron jumps across 
I Landau level spacings with simultaneous absorption of 
m photons of frequency w\ and emission of ri2 photons of 



5 



frequency u>2, or simultaneous emission of n\ photons of 
frequency w\ and absorption of n 2 photons of frequency 
u>2 (thus the symbol 77.2072 — n\U)\ : I has the same mean- 
ing as n\tj)\ — ri2o;2 : V). The symbol n\u>\ : I indicates a 
process during which an electron jumps across I Landau 
level spacings with the assistance of n\ real (emission or 
absorption) 07i-photons and virtual 072-photons. So does 
the symbol 772072 '■ I- We find that the processes repre- 
sented by n\u>\ + 77,2^2 : I contribute, in the R xx -vevsus- 
o7 c /o7 curve, a structure consisting of a minimum and a 
maximum on both sides of O7 c /o7i = n^[l + 712072/771071)// 
or lu c /uj2 = 1^2(1 + niLUi/n 2 L02)/l- And those by 7ii07i — 
U2UJ2 '■ I (assume n\UJi > 772^2) contribute a minimum- 
maximum pair around lo c /lo\ = m(l — 772072/771071)// or 

W c /w 2 = 77 2 (771071 /772W2 — I)/'- 

From the structure of R xx curves with increasing ra- 
diation strength in Fig. 1 we can clearly identify sev- 
eral peak-valley pairs characteristic of the bichromatic 
irradiation. The strongest peak-valley pair is around 
w\/u) c ~ 2 or 072 /o7 c ~ 3, which can be referred to the 
joint contribution from single and multiple photon pro- 
cesses with mono-07i, mono-W2 and mixing u>i and IV2, 
such as o;2 — 07 : 1, 07i : 2, 072 : 3, 2o7i : 4, u>i + u>2 : 5, • • • 
. The other peak-valley pairs which are clearly identified 
include that around 07i/o7 c ~ 0.4 referred to u)\ + 072 : 1, 

• • • ; that around o7i/o7 c ~ 0.57 referred to 2o7i + 072 : 2, 

• • • ; that around O7i/o7 c ~ 0.8 referred to 071+072 : 2, 

• • • ; that around o7i/o7 c ~ 1.14 referred to 2o7i + 072 : 4, 

• ■ • ; that around o7i/o7 c ~ 2.84 referred to 2o7i +072 : 10, 
■ ■ • . They are indicated in the figure. The minima of all 
these peak-valley pairs can drop down to negative when 
increasing the strengths of the radiation field, indicating 
the the possible locations of the bichromatic ZRSs. 



To compare the photoresponse under bichromatic ra- 
diation with those under monochromatic radiation we 
plot in Fig. 2 the magnetoresistivity R xx versus the in- 
verse magnetic field for the above system under simul- 
taneous irradiation of two microwaves having frequen- 
cies 07i/27r = 31 GHz and o72/27r = 47 GHz with in- 
cident amplitudes En — E\2 — 3 and 3.5V/cm, to- 
gether with those subjected to single microwave radia- 
tion of frequency 07i/27r = 31 GHz with incident ampli- 
tude E[ = 3 V/cm, or of frequency uji/2-k — 31 GHz with 
incident amplitude E\ = 3 V/cm. 

This figure is redrawn in Fig. 3 showing R xx as a func- 
tion of the magnetic field to give a clearer view in the 
higher magnetic field side. It indicates a possible bichro- 
matic ZRS around O7 c /o7i ~ 2.2 or 07 c /o72 ~ 1.5 arising 
from R xx droping to negative at the minimum of the 
valley-peak pair around O7 c /o7i ~ 2.5 associated with the 
mixing biphoton process 071+072 : 1. This result is in 
agreement with the recent experimental observation on 
the possible precursor of bichromatic ZRSii^ 

Note added: After the acceptance of this paper by 
Phys. Rev. B, we read a further report of the bichro- 
matic microwave photoresistance measurement^ The 
experimental peak/ valley positions exhibit a good agree- 
ment with the present theoretical results in a wide mag- 
netic field range. Our predicted peaks around O7 c /o7i = 
0.77,0.52,0.37, and valley at 07J071 = 0.33 (Fig. 3), 
quite accurately (to = 0.068 m c ) reproduce the observed 
peaks around B w 570, 380, 275 G and the valley around 
B « 250 G. 

This work was supported by Projects of the National 
Science Foundation of China and the Shanghai Municipal 
Commission of Science and Technology. 
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